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This paper presents a structural dynamic modeling of the Long Duration Exposure 
Facility (LDEF) , which is a Space Shuttle payload of passive scientific experi- 
ments contained in trays mounted on a large cylindrically shaped structure. 

Special detailed finite element modeling, using the SPAR system of computer pro- 
grams, was required to obtain good agreement between analytical and test vibration 
modes. The scientific experiment trays contributed significantly to overall LDEF 
stiffnesses, and these contributions were realistically represented for each tray 
by the stiffness matrix of an equivalent orthotropic panel in the overall LDEF 
model. Orthotropic stiffnesses for this equivalent panel were obtained from 
finely detailed statically loaded SPAR models which included stiffness coupling 
and allowed for partial relative sliding of the tray clamping attachments. Sen- 
sitivity to LDEF joint boundary conditions was determined, and static test data 
proved valuable in assessing modeling of local end fittings. 
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SUMMARY 


The purpose of this paper is to present the recent application of structural 
modeling of the Long Duration Exposure Facility (LDEF) utilizing the SPAR 
system of computer programs for vibration analysis. 

The technical areas of interest consist of five parts: 

1. Development of the LDEF finite element model. 

2. Derivation of tray effective panel stiffness matrix using finite element 
tray models. 

3. Assessment of attachment conditions and end fitting flexibility by com- 
paring SPAR with test static displacements. 

4. SPAR Grouping. 

5. Derivation of the LDEF frequencies and mode shapes and comparing 
them with test. 

Results indicate that the SPAR LDEF model gave good agreement between 
analytical and measured data. This was true for the frequencies of the LDEF 
model without trays and the LDEF model with trays. The static test data is 
valuable in assessing end fitting boundary conditions to be used in the finite 
element model. 
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INTRODUCTION 


The Long Duration Exposure Facility (LDEF) vehicle will be launched 
in the space shuttle and placed in earth orbit with 86 trays, most of them con- 
taining passive experiments. After a period of time exposed to the space 
environment, LDEF will be retrieved by a second space. shuttle flight, and 
its many experiments will be evaluated. 

Because LDEF must interface with the space shuttle, it is important 
that the LDEF finite element dynamic model be an accurate representation of 
the actual structure. NASA Langley performed dynamic model survey testing 
and static testing of LDEF with and without the trays being installed. The 
SPAR finite element computer programs were used to model LDEF and the 
static test data was used for correlation and model improvement. SPAR has 
the advertized capability to perform a vibration analysis of a system with 
thousands of degrees-of-freedom (DOF) but the analyst must check eveiy 
input to maintain the accuracy desired. The philosophy used in this LDEF-SPAR 
model was to keep the model as simple as possible but maintain the required 
accuracy and utilize a thorough checking of inputs. Each experiment tray was 
modeled as a single effective fully coupled six-by-six stiffness matrix panel. 

The stiffness matrix was obtained from static loadings on SPAR tray models. 
One of these tray models has 4356 DOF compared with 1416 DOF for the 
entire LDEF structure. Therefor, with 86 trays to model, the effective panel 
approach becomes most attractive. 

The information in the body of this report presents the LDEF model 
construction, static test correlation for determining the joint boundary con- 
ditions and end fitting flexibility, how the tray effective panels were derived, 
and finally the comparison of frequencies and mode shapes with test data. 

The information in the appendices consists of the panel and beam prop- 
erties, the mesh size criterion, the tray moments of inertia, and the SPAR 
eigenvalue extraction experience. 
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SYMBOLS 


General Notations 
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Tray effective panel stiffness term - force in i direction 
due to a unit strain in j direction 

distance between two panels 

tray panel offset due to longitudinal extensional-bending 
coupling 

tray panel offset due to lateral extensional-bending coupling 

tray panel offset due to shear-twisting coupling 

single panel bending stiffness 

effective panel bending stiffness 

extensional strain 


shear strain 


Young's Modulus 
force 


free boundary condition, applied load column matrix 
Shear modulus 

thickness, height of circular arch at mid point 
total tray plus ballast mass moment of inertia 
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tray ballast mass moment of inertia 

tray mass moment of inertia 

rotational strain 

rotational strain in x direction 

rotational strain in y direction 

twisting strain 

stiffness , stiffness matrix 

stiffness of rigid element 

length of effective tray panel 

moment 

bending stress resultant 

reaction moment at grid point i in the tray coordinate 
direction j 

extensional stress resultant 
applied static test load 

reaction force at grid point i in the tray coordinate 
direction j 

rotation 

static test reaction load 
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simple support boundary condition 


thickness 


tray panel longitudinal effective thickness 


tray panel lateral effective thickness 


tray panel shear effective thickness 


extensional displacement 


displacement column matrix 


weight, deflection normal to surface 


weight of tray ballast 


width of effective tray panel 


tray coordinate system 


tray ballast mass moment arm 


tray mass moment arm 


percent displacement error by adding a rigid link in 
series with a flexible element 


deflection 


rotation displacement 
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Poisson's Ratio 




(0 


stress in i direction 


angle of rotation between tray coordinate system, x,y, z, 
and the LDEF SPAR cylindrical coordinate system, 
r 0 z' 


SPAR Notations 

ANG angle section 

BA general beam elements 

BB beam 6x6 stiffness matrix 

E21 general beam elements such as angles, tees, zees, 

tubes, etc. 

E25 zero-length element used to elastically join two coincident 

grid points 

E33 triangular membrane plus bending element 

E43 quadrilateral membrane plus bending element 

BRL beam rigid link offset 

CHN channel section 

GIVN beam with extensional and bending stiffness specified 

SA shell section properties 
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TUBE 


tube section 


TEE 


tee section 


WFL 


I section 


Experiment Tray Notation 

DP 7. 62 cm (3 in) side Debris experiment tray 

S 7. 62 cm (3 in) side shallow tray 


N 


15.24 cm (6 in) side nominal tray 


D 


30. 48 cm (12 in) side deep tray 


SR 


7.62 cm (3 in) end shallow rectangular tray 


NR 


15.24 cm (6 in) end nominal rectangular tray 


dr 


30. 48 cm (12 in) end deep rectangular tray 


SC 


7.62 cm (3 in) end shallow square (comer) tray 


NC 


15. 24 cm (6 in) end nominal square (comer) tray 


DC 


30. 48 cm (12 in) end deep square (comer) tray 


7.62 cm (3 in) rectangular end shallow tray stiffness 


15.24 cm (6 in) rectangular end nominal tray stiffness 


^12 


30.48 cm (12 in) rectangular end deep tray stiffness 



LDEF STRUCTURAL TEST CONFIGURATION 


The Long Duration Exposure Facility (LDEF) structure is primarily a 
bolted assembly of 6061 - T6 aluminum extrusions. The LDEF structure 
schematic is shown in Figure 1 indicating the principal dimensions and the 
design features for the 12 sided regular polygon. The polygon is 4. 267 meters 
(14 feet) across the points and 9. 144 meters (30 feet) long. The center and 
end bulkheads anchor the longerons and provide reaction points for the 8 
diagonal tubes. The intercostals connected to the longerons divide the per- 
iphery into 72 equally sized rectangles. The center and end bulkheads are 
welded, the other joints are bolted. The center bulkhead has 2 trunnion 
fittings and a keel fitting that support the LDEF in the Shuttle bay. The 
forward bulkhead carries a trunnion fitting which acts thru a cross beam 
(Figure 3) to provide 2 other attachment points to the Shuttle. Figure 2 shows 
the aft end of the LDEF dynamic test configuration with the simulated trays 
installed. Figures 3and 4 show the forward end and the right side respectively 
of the LDEF structure without trays in the transportation vehicle. 

The longeron-intercostal joint design is identified in Figure 5 by the 
triangular arrows. These are pin-clevis joints at each end of the intercostals 
as shown in Figures 12 and 13. The clevis bolts are in the LDEF longitudinal 
direction and can be considered as pinned joints in a finite element model. How- 
ever, the clevis tolerance is such that some relative rotational motion could occur 
about the 2 mutually perpendicular axes normal to the bolt axis. Therefor, the 
analysis considered both pinned and ball joints at these connections. 

The eight diagonal tubes can best be seen in Figure 4 with the detail of 
the end joint at the center bulkhead shown in Figure 6. These clevis type pin 
joints are in reality ball joints because of their internal design. Two ball joints 
connect to a common flat plate type fitting extending from the center bulkhead. 

The original SPAR model considered this flat plate fitting as a rigid link while 
a later refinement included its flexibility. 

The LDEF side trays are riveted 6061 - T6 aluminum alloy angle and 
sheet of either 0. 160 cm (.063 in) or 0. 3175 cm (0. 125 in) thickness with 
0. 3175 cm (. 125 in) Tori section bottom cross beam members. Both the T and 
I sections can be seen in Figure 7. The side tray experiment envelope is approx- 
imately 1. 219 m long and 0. 914 m wide with depths of 7. 62, 15. 24 and 30. 48 cm (3,6 
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and 12 in). The end tray envelope is approximately . 762 m square with depths 
of 7.62, 15.24 and 30.48 cm (3, 6 and 12 in). Figure 7 shows a 7.62 cm (3 in) 
side shallow tray mounted on a longeron-intercostal frame. The base plate of the 
experiments are bolted to the bottom framework at some of the many nut plate 
locations provided. For the LDEF dynamic configuration, all of the experi- 
ments except the Debris tray were represented by a single 0.457 cm (0. 18 in) 
aluminum base plate with 1, 2 or 3 stacks of aluminum ballast plates to repre- 
sent the correct experiment mass. Each ballast stack was bolted to the base 
plate at the four comers of the stack. These bolts can be seen in Figure 2 and 6. 
The Debris tray consists of two 0.457cm (0. 18 in) aluminum plates with no ballast 
plates, with each plate covering one half of the tray bottom. Each Debris plate 
was bolted to the tray bottom frame with 15 bolts. All of the other trays with a 
single base plate were bolted to the tray bottom with 6 bolts, 4 at the comers and 
2 to the center framework. Figures 7 and 8 show the tray stmcture interface for 
a side tray. There are no bolts passing through the tray lip. The tray clamp 
design allows the tray to expand in the plane of the tray, due to thermally induced 
loads,by overcoming the clamp friction load. For the dynamic analysis the four 
comer clamps were considered to have no relative motion between the tray lip 
and the LDEF stmcture. The four side clamps were assumed to allow relative 
rotation about an axis normal to the tray lip. 
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Rom Position; 


Bay Positions 



Note; 

1. Location for a Magnetic Viscous Damper. 

2. All openings on the LDEF are covered by either a tray or plate. 

3. Longerons - 17. 78 cm (7 in) Aluminum I-Section, End Frames - 15. 24 cm (6 in) Aluminum I-Section, 

Diagonals - 15.24 cm (6 in) Aluminum Tube. Intercostals are Aluminum T-Sections, one half 
of a 17.78 cm (7 in) I-Section. 

4. Auxiliary Trunnion for Ground Handling Only. Removed before flight. 


Figure 1. LDEF Structure Schematic 



















Figure 8. Tray Structure Interface 


SPAR LDEF MODEL 


The SPAR computer programs of Reference 1 were used for the LDEF 
finite element modeling and for the subsequent static and dynamic analysis. 

Figure 9 shows a computer print out of the undeformed shape of the SPAR model. 
The bottom, side and end views are shown. The model has 236 grid points with 
1416 degrees of freedom. There are 381 beam elements, 16 three node elements 
and 92 four node elements. The circles in Figure 9 represent the 16 ball joints 
at the ends of the 8 diagonals. The model shown is the final configuration that 
includes the diagonal end flexibility of the flat plate fittings between the center 
bulkhead and the ball joints. The original model had rigid links for these flat 
plate fittings. The triangle represents a pin joint that allows the cross beam to 
rotate about the 3 axis relative to LDEF. The dots show the longeron- to-inter- 
costal joint locations that were treated as pin or ball joints in the analysis. In 
most cases each dot represents two pin joints as shown in Figures 12 and 13 so 
that the 12 grid points per station results in 24 pin joints per station. However, 
at the forward and aft bulkheads.the top and bottom intercostals are of different 
construction and do not have clevis pin joints as shown in Figure 5. There are 
still 12 grid point locations but only 20 pin or ball joints per station making a total 
of 136 of these joints in the LDEF model. There are no pin or ball joints at the 
center bulkhead. 

SPAR has the capability of subdividing the model into groups of substructure 
elements. This grouping has the advantage of visually checking each group plot for 
modeling errors and obtaining the weights of each group for comparison with 
actual weights for weight distribution control. The SPAR model was divided 
into the substructure groups shown in Table 1. 
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Figure 9. LDEF SPAR Model 
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TABLE 1. SUBSTRUCTURE GROUPS 


Modeling Sequence 

SPAR Group 

1 

Longerons 

E21-1 

2 

Intercostals 

E21-2 

3 

Center BHD 

E21-3 

4 

Forward BHD 

E21-4 



E33-1 

5 

Aft BHD 

E21-5 



E33-2 

6 

Diagonals 

E21-6 



E25-2 



E21-9 

7 

Cross Beam 

E21-7 



E25-1 

8 

lECM Shelf 

E21-8 



E43-1 

9 

Exp. Trays 



Debris Side 

E43-2 


Shallow Side 

E43-3 


Nominal Side 

E43-4 


Deep Side 

E43-5 


End- Rectangular 

E43-6 


End-Comer 

E43-7 


Element Type 

No. of Ele. 


Beam 

96 

10 

Beam 

72 

11 

Beam 

30 

14 

Beam 

40 

15 

Triangular Plate 

8 


Beam 

40 

16 

Triangular Plate 

8 


Beam (Tube) 

32 

17 

Zero length beam 

16 

9 

(ball joints) 

Beam (end attach. ) 

20 

17 

Beam 

13 

18 

Zero length beam 

1 

9 

(pin joint) 

Beam 

21 

19 

Quadrilateral Plate 

6 



Quadrilateral Plate 

22 

Quadrilateral Plate 

26 

Quadrilateral Plate 

14 

Quadrilateral Plate 

10 

Quadrilateral Plate 

6 

Quadrilateral Plate 

8 
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Longerons and Intercostals 

The first substructure group modeled was the longeron group. The 
longerons consist of 12 forward and 12 aft I beams attached to the center bulk- 
head. Each longeron, at location A through L in Figure 10, was divided into four 
beam elements as shown at the intercostal attachment points resvilting in a total 
of 96 beam elements. Figure 11 identifies all of the longeron grid point numbers 
and the SPAR section property used. The second substructure group modeled 
consisted of the 72 intercostals tabulated in Figure 11. Figures 12 and 13 show 
the actual longeron to intercostal attachment and the SPAR model configuration. 

The longeron beam elements were connected to grid points located on the lor^eron 
neutral axis. The intersection of the neutral axes of two adjacent intercostals 
fell outboard of the longeron neutral axis as shown in Figure 12. The intercostal 
beam elements were modeled with beam rigid link offsets (BRL) to locate the 
intercostal beam elements at their neutral axes. To avoid adding two additional 
grid points attached by a zero length beam elements for each pin joint, the inter- 
costal beam elements bending stiffness was set equal to zero to simulate the pin 
ended condition. For a ball-joint end conditions both of the bending and the tor- 
sional stiffnesses of the beam elements were set equal to zero. The SPAR 
section property BA16 was used for the ball-joint end condition and BA18 
was used for the pin-joint end condition. Both of these section properties 
have some of the tee section properties (BA2) set equal to zero as shown in 
Appendix A (Table 9). 
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( 2 ) 


Center Bhd. 



A Longeron Grid Point Locations 

* ^ I (at Long. N. A. , Eadius = 2.0447) 

• B 

» C 


'd 

View Looking Aft 

I 

E 

( 

F 


( 2 ) 

L. 

K f 

J* 

( 1 ) + 

I* 

t 

H 

I 

G 


Location Coordinates 



(1) 

(2) 

A 

-. 529206 

1. 975028 

B 

-1.445821 

1.445821 

C 

-1.975028 

. 529206 

D 

-1.975028 

. 529206 

E 

-1.445821 

-1.445821 

F 

-. 529206 

-1.975028 

G 

. 529206 

-1.975028 

H 

1.445821 

-1.445821 

I 

1.975028 

-. 529206 

J 

1. 975028 

. 529206 

K 

1.445821 

1.445821 

L 

. 529206 

1. 975028 


Figure 10. Longeron Grid Point Coordinates (All Dimensions in Meters) 
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Fwd Bhd. 

-4.4313 

Intercos. 

-4. 3498 

Intercos. 

-2.9401 

Intercos. 

-1. 5304 

Center Bhd. 

0.0 

Intercos. 

1. 5304 

Intercos. 

2.9401 

Intercos. 

4. 3498 

Aft. Bhd. 

4.4313 


A 

B 

C 

Longeron Location (See Figure 10) 

D E F G H I J 

K 

L 

Grid Point Numbers 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

113 

114 

115 

116 

117 

118 

119 

120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

i 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 


No. of 
Elements 

B'Al 96 17.78 cm (7. in) I beam - all lor^erons 

BAl 4 top and bottom intercostals at end stations 

BA16 68 Ball Joint^ 

or / 8.89 cm (3. 5 in) T section (BA2) - all other intercostals 

BA18 68 Pin Joint 


SPAR Section 
Property 
(Table 9) 


Figure 11. Longerons and Intercostals 
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Actual Configuration 


Longeron 



Intercostal Beam 
Element 



Common Pin or Ball Joint 


Intercostal Beam Element 
(at Neutral Axis) 


Intercostal Beam Rigid Links 
(BRL) 

4 

Grid Point 


(At Longeron Neutral Axis) 


Figure 12. Longeron-Intercostal Attachments at Stations ± 4. 4313m 
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I 


Actual Configuration 



SPAR configuration similar to Figure 12. 


Figure 13. Longeron Intercostal Attachments at Stations ± 1. 5304m and ± 2. 9401m 
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Bulkheads 


The third substructural group modeled was the center bulkhead made up 
of 30 beam elements as shown in Figure 14. There are 14 vertical and lateral 
I beam members with the outer 14 ring elements representing a box section made 
up of two channel sections with inner and outer webs attached. Grid points 71 
through 82 are the longeron grid points. Grid points 83 and 88 are the locations 
of the main support trunnion fittings and grid points 84 and 87 the attachment 
points for the 8 diagonals. Grid points 99 and 100 are the attachment points for 
each side of the lECM shelf. All beams are located at their neutral axes and 
beam rigid link offsets were used to join the beam element ends. The 
fourth and the fifth substructure groups are the identical forward and aft bulk- 
heads shown in Figures 15 and 16. They each have 40 beam elements and 8 
triangular plate elements. Grid points 10 through 21 and 149 through 160 are 
the longeron grid points. Grid point 28 connects to the cross beam. All of the 
vertical and horizontal beam elements represent I beam sections except the four 
connected to the center grid points. These four beams are box sections made up 
of two channel sections opening outboard with plates on the forward and aft sur- 
faces. The eight triangle plate elements represent two plates 15.24 cm (6 in) 
apart. These two plates were input as a single effective panel stiffness matrix by 
using the FORMAT = 2 option in shell section properties (SA) inputs. The 
derivation of these properties is shown in Appendix A. The two plates are 
assumed to bend as a single panel and have an extentional stiffness equal to the 
sum of the two plates. The effective panel bending stiffness is: 


D 


Et 


2 

(1-v) 


( i d^ + dt + |- t^) 


where 


t = . 3556 cm (. 14 in) 
d = 15. 24 cm (6. in) 


( 1 ) 
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SPAR 

Saction 


(Table 9) 


BAl 

BA3 

BA4 

BAS 


No. of ♦ 

Elements 

5 - 17. 78 cm (7 in) I beam 
3 - 30. 48 cm (12 in) I beam 

6 - 24. 13 cm (9. 5 in) I beam (average section) 
16 - Built up Box Section 


Figure 14. Center Bulkhead at Station 0. 0 
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5062m 



(36 places) 


SPAR 

Section 


No. of 
Elements 


Beam Rigid links (BRL) 
(8 places) 


36 - 15. 25 cm (6 in) I beams 
4 - 35. 56 X 25. 4 cm (14. x 10. in) built up box beams 
8 - Triai^ular panel elements 


Figure 15. Forward Bvilkhead at Station - 4. 4313m 









Diagonals 


The sixth substructure group modeled was that of the eight diagonal tubes with 
the end ball joints and end attachment fittings shown in Figure 17. The circles 
represent the 16 ball joints. Flexible beam elements connect the ball joints to 
grid points 84 and 87 on the center bulkhead and at eight grid points on the 
longerons near the end bulkheads. The original model had rigid links from 
the ball joints to these end grid points. The ball joints were modeled by locating 
two grid points at the ball joint location and connecting them with a zero length 
element (E25) using a beam six-by-six intrinsic stiffness matrix input (BB) in 
SPAR. The stiffness matrix, K , is defined as 

F = KU (2) 



In the above f, m , u and r equal applied force, applied moment, displacement 
and rotation at the origin of the beam element in the global coordinate system. 
For the diagonal ball joints, the K matrix takes the form 



30 



The diagonal term is set equal to zero for each unrestrained degree of freedom 
desired. The diagonal terms not set equal to zero represent the'rigid spring*(Kp 
that is in series with the flexible element. The value must be large enough to 
limit the displacement error but not so large as to result in computer computa- 
tional errors. The'rigid link stiffness',' , required to obtain a given displace- 
ment error ( a ) was obtained from: 


lOOK 



where 


^4) 


a 


The percent displacement error of a"rigid"element (Kp 
in series with a flexible element (K) compared with the 
displacement of K by itself. 


For the diagonal elements, the extensional stiffness, K , of the upper diagonal 

6 9 

was 48.6906 MN/m (.27803 x 10 Ib/in). The value of 700508 MN/m (4.0 x 10 

Ib/in) was used for indicating a diagonal extensional error of 0. 007 percent. 

This value of K was used for K , K and K in equation (3). 

X* XX (JO 

Each tube was divided into four equal beam elements. The mesh size 
criterion is discussed in Appendix B. The criterion used states that the natural 
frequency of a simple supported beam element must be greater than 50 Hz. If 
this criterion was not met an intermediate grid point was added. For the diagonals, 
with the length between ball joints used, the natural frequency was 20. 5 Hz and with 
one intermediate grid point at the center the frequency was 82 Hz. The additional 
quarter point locations were added because of the actual ball-joint boundary con- 
ditions at the diagonal ends ; the criterion is actually intended for structures 
with built in ends. 
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Lateral Locations (1) GP 193-212 at ± 1. 8796m 

GP 185,186,191,192,213,214,219 & 220 at ± .63391m 
GP 187-190,215-218 at ± .61468m 


( 3 ) 


( 3 ) 


SPAR 
Section 
Properties 
(Table 9) 

BA8 

BA19 

BA20 


No. of 
Elements 


32 - 15. 24 cm (6 in) tube (all main diagonal elements) 

4 - flat plate 2. 54 x 18. 288 cm (1. 0x7.2 in) with 18. 542 cm (7. 3 in) 
BRL at grid point 84 and 87. 

16 - flat plate 1. 905 x 17. 78 cm (. 75 x 7. 0 in) connecting ball joints 
to grid points 195,200,205,210,35,40,41,46,137,142,143 & 148 


Figure 17. Diagonals 
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Cross Beam 


The cross beam shown in Figure 18 was the seventh substructure group 
modeled. This group consists of 13 beam elements and one zero length beam 
element, E25 , between grid points 9and 184 to represent the pin joint. The 
spar input stiffness matrix , BB, was of the form of equation (5) for a pin joint 
about the 3 axis. All values on the diagonal of the stiffness matrix were given 
a value of 700508 MN/m or a value at least two orders of magnitude larger than 
the bending stiffness of the spindle pin (3 elements connecting grid points 184 
to 28). The error in the bending deflection of a simple cantelever beam in 
series with a rigid link using equation (4) was less than one tenth of one percent. 
The four lateral tubes have relatively stiff end fittings. The fittings were made 
rigid by using beam rigid link offsets so that the flexible tube lengths were 
correct. The upright beam elements connecting grid points 3 to 9 and 7 to 9 
are of a complex I beam configuration seen in Figure 3. The torsional stiffness 
due to differential bending was added to the basic torsional stiffness of these 

4 

beam elements. The basic section had a torsional stiffness of 105.7 cm com- 

4 

pared with the differential-bending torsional stiffness of 1111. cm . This 
stiffness controls the frequency of the first flexible mode. 

The SPAR input stif&iess matrix, BB , was of the form 


K 


K 

0 

0 

0 

0 

0 


11 


0 

K 

0 

0 

0 

0 


22 


0 

0 

K 

0 

0 

0 
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0 

0 

0 

K 

0 

0 


44 


55 


0 

0 

0 

0 

0 

0 


for a pin joint about the 3 axis. 


( 5 ) 
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= beam rigid link (BRL) 


(2) 



Pin Joint 
About 
3 Axis 


SPAR 

Section 


Properties 


(Table 9) 


No. of 
Elements 



BA9 8 - lateral members are 10. 16 cm (4 in) dia. tubes 

BAIO 2 - built up section, torsion stiffness includes differential bending 

BAll 1 - complex flexible section at center of spindle pin 

BA15 2 - solid circular shaft 


Figure 18. Cross Beam at Station - 4.4313m 
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lECM Shelf 


The eighth substructure group was the lECM shelf structure consisting 
of 21 beam elements and 6 quadrilateral plate elements as shown in Figure 19. 
Grid points 85, 86, 99 and 100 are located on the center bulkhead. The 
5. 08 cm (2. in) diameter diagonal tubes support the lECM shelf off of the 
center bulkhead. The outer frame and the support beams at stations . 2794 m 
and . 9906 m are made of channel sections. The other panel stiffeners at 
stations . 5364 m, .7396 m and 1.2002 m are angle beam elements. There 
are six . 160 cm (. 063 in) thick quadrilateral plate elements. 
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SA2 



SPAR 

Section No. of 

Properties Elements 
(Table 9) 


BA12 

BA13 

BAH 

SA2 


16 - 7.62 X 4.445 cm (3. x 1.75 in) channel 
3 - 2. 54 X 2. 2352 cm (1. x . 88 in) angle 
2 - 5. 08 cm (2. in) dia. tube 
6 - . 160 cm (. 063 in) plate 


Figure 19. lECM Shelf 
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Experiment Trays 


The ninth and final substructure group included the 86 experiment trays that 
were modeled as equivalent orthotropic panel elements as depicted in Figure 20. 

A static SPAR tray model was used to obtain the six by six panel stiffiiess 
matrix for each type of tray. Figure 21 shows the tray locations. 


Figure 21 


Types of Tray 

Tray Depth 

Code 

Side Trays 

Debris experiment 

7.62 cm (3 in) 

DP 

Shallow 

7.62 cm (3 in) 

S 

Nominal 

15.24 cm (6 in) 

N 

Deep 

30.48 cm (12 in) 

D 

End Trays 

Shallow rectangular 

7. 62 cm (3 in) 

SR 

Nominal rectangular 

15. 24 cm (6 in) 

NR 

Deep rectangular 

30.48 cm (12 in) 

DR 

Shallow square (comer) 

7. 62 cm (3 in) 

SC 

Nominal square (comer) 

15. 24 cm (6 in) 

NC 

Deep square (comer) 

30.48 cm (12 in) 

DC 


As shown in Figure 20 there is coupling between extensional and bending terms as 
well as between the shear and twisting terms. Figure 22 shows the extentional- 
bending coupling of a shallow side tray. Figure 23 shows the shear-twisting 
coupling of the base plate of a deep side tray. In Figure 21 each rectangle 
represents a tray with the type of tray in the upper left hand comer and the 
total tray weight in the lower right hand comer. The bay location and the row 
numbers are also shown in Figure 21. The tray model shown in Figure 20 is the 
Debris experiment tray. All of the four side tray types were modeled and the 
nominal and deep rectangular end trays were modeled. 
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Longeron 





4 101.60 36.29 


; 24.95 100.70 

4 - 

!DP 1 


41.28^ 

j 101.60 


DP 

100.70 

24. 9J 


24.95 


101.60 


101.60 


63.05 


Fwd End 


102. 06 


102. 06 



E 

F 

S 

44.45 

S 

44.45 




S 

102.06 

DP 

24.95 

N 

101.60 

DP 

24.95 

S 

102.06 

"dp 

24.95 

DP 

51.71 




DP 

24.95 


39.92 

N 

41.73 

"dp" 

24.95 

s 

72.12 



102.06 


163.75 


24.95 


Aft End 



Figure 21. Tray Locations and Total Weight (kg) (See Figure 1) 













































Figure 22. Static Deformation of Tray Under Unit Longitudinal Strain 
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Figure 23. Shear- Twisting Coupling of the Deep Side Tray Base Plate 

(Shear Strain Equal to Unity) 



TRAY EFFECTIVE PANEL STIFFNESS PROPERTIES 


Stiffness Matrix From SPAR Tray Models 


The flanges of the side trays are clamped to the structure at 8 locations 
as shown in Figure 8. Rigid links were used to connect these 8 clamping points 
on the tray lip to the plane of the four comer longeron grid points and to four 
points mid way between these grid points. 
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Tray grid points 1 , 3, 4 and 6 are at the actual location of the LDEF longeron 
grid points. Tray grid points 2, 5, 7 and 8 are located midway between the 
longeron grid points. 

The SPAR panel stiffness matrix is of the form 
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Si 

S2 

0 

^24 

^25 
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N 

xy 


0 
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•=33 

0 

0 

Se 


e 

xy 

M 

X 


Si 

S 2 

0 

S4 

Ss 

0 


k 

X 

M 

y 


Si 

^52 
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^54 

^55 

0 
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y 

M 

^ xy> 
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0 

S3 
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0 

S6_ 


k 

V xy^ 


42 


where 


3 w 


3 w 


k = 

X 


3 w 

2 ’ 

3 X 


3y 


k 

xy 


2 

3x 3y 


Each column of this stifftiess matrix was obtained by setting one of the strains 
equal to unity and all others equal to zero. The six resulting conditions are 
the longitudinal strain equals unity ( e =1. 0), the lateral strain equals unity 

X 

( e^ = 1. 0) , the shear strain equals unity ( e = 1. 0), the longitudinal rotation 

equals unity ( k =1. 0), the lateral rotation equals unity ( k =1. 0), and the 

X y 

twisting rotation equals unity ( k^ = 1. 0). Table 2 shows the actual displace- 
ments imposed on the tray model at the ends of the rigid links for each of the 
six displacement conditions. The force reactions output at grid points 1 through 
8 were used to determine the stiffness coefficients in equation (6). The stress 
resultants on the left side of equation (6) can be expressed as functions of these 
reaction forces. 
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(7) 


where 

P.. Reaction Force 
1 ] 

M.. Reaction Moment 

13 

Subscript i is the grid point number 
Subscript j is the tray coordinate direction 
L Length of effective panel 
W Width of effective panel 

From equation (6) and the boundary conditions in Table 2 the values of the stiffness 
matrix are calculated by inserting the reaction forces and moments into equation (7). 
The first column of the stiffness matrix is obtained by setting e = 1. 0 in 
equation (6) and all others equal to zero. 



( 8 ) 
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1.0 


1.0 


for e 


(n) = 


xy 


for k = 

X 




0 


^14 

0 


S4 



0 

33 

0 


c.. 

44 

0 


C 

c 


0 

63 







The design of the tray-to-structure clamp allowed in plane relative sliding motion 
between the tray flange and the mounting structure. The four side mid point 
clamps were allowed some relative motion by letting 0 ^ be free to rotate at 
grid point 2, 5, 7 and 8. This is shown in Table 2 for the shear (e^ = 1. 0) 
and the twisting = 1.0) boundary conditions. By allowing this flange motion 
the shear and twisting stiffnesses are reduced to 54 and 58 percent of the values 
for the original boundary conditions, as shown in Appendix A. 

To demonstrate the procedure, the numerical results from the evaluation 
of the Debris tray subjected to a longitudinal strain equal to unity ( 
be shown in the following discussion. The reaction forces and moments from the 
SPAR static run are shown below. 


REACTION FORCES (MN) REACTION MOMENTS (Mm-N) 
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With a panel length (L) of 1. 5304 m (60.25 in) and a panel width (W) of 1.0584 m 
(41.6698 in) the stress resultants from equations (7) and (8) result in the first 
column of the stiffiiess matrix 
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X 

o 

11 

. 3604 MN/cm 
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M 
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II 
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4. 034 MN 

M 

y 

o 
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II 

-.4675 MN 

M 

xy 

Si = 

0.0 


All of the other five static boundary displacements are imposed on the same tray 
model to obtain the six-by-six stiffiiess matrix. 


Original Debris Tray Stiffness Matrix (MN/cm, MN & Mm - N) 
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Because some of the off diagonal terms are not equal, the average values were 
taken resulting in the final stiffness matrix for the Debris tray of 
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The final stiffness values used in all the equivalent panels are tabulated in 
Appendix A. 

The units for the stiffness matrix are partitioned as 
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TABLE 2. DISPLACEMEMT BOEINDARY CONDITIONS 
FOR TRAY STATIC MODELS 
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Extrapolation of Rectangular End Tray Stiffness 

Because the end bulkhead structure is of a rigid welded I beam configur- 
ation, the end tray panel stiffness will have a minor effect on the LDEF dynamic 
response. For this reason only two of the end trays were modeled, the nominal 
and deep rectangular trays, and the stiffness for the other trays were extrapolated 
from this data. The valueof each stiffness term of the shallow tray rectangular 

matrix ( C ) was extrapolated linearly from the nominal ( C. ) and deep ( C,_ ) 
o o 12 

rectangular stiffness terms: 

C 3 ' C6-<S2-06> = 2 O 3 - (9) 
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Tray Panel Effective Properties 

The construction and clamping of the rectangular and square end trays 
were very similar so the following set of equations were used to obtain the 
square trays stiffness matrix knowing that the x and y direction properties 
must be equal. The effective tray panels can be represented as orthotropic 
panels that have membrane-bending coupling and twisting - shear coupling be- 
cause of rigid linlc offsets. 

, oirthotropic panel 



The orthotropic stiffness matrix has the form: 
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Because of the offset ( d ) there are coupled loadings in the xz plane. 
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Using these coupling equations, the left side of equation (6) can be rewritten. 
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By equating the first row equation with the fourth row 


11 

o 

to 

0 

^14 

®15 

0 

21 

^2 

0 

^4 

^5 

0 


0 

^33 

0 

0 


41 

^=42 

0 

^44 

‘^45 

0 

51 

°52 

0 

^54 

^55 

0 


0 

^63 

0 

0 

^66 


X 

e 

y 

e 

xy 

k / 

X 

k 

y 

k 

xy 




d = 

X 


C.^e + C._e + C..k + C._k 
41 X 42 y 44 x 45 y 

+ C,„e + C, .k + C, ^k 
11 X 12 y 14 X 15 y 


By equating the second row and fifth row equation 


d = 

y 


c e + C__e + C_.k + C__k 
51 X 52 y 54 x 55 y 

C_,e + C_„e + C-.k + C„^k 
21 X 22 y 24 X 25 y 


By equating the third row and sixth row equation 
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The three load coupling offsets d , d and d change with changes in the 
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boundary strain conditions. For extrapolating purposes these equations were 



reduced to the following by using the boundary conditions : 
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Using the orthotropic equations for composites from Reference 3, the values of 
the Poisson's Ratios and the extensional and shear effective thicknesses for each 
type of panel can be derived. 
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From equations(14) and (15) the orthotropic properties and the offsets of each tray 
type are tabulated in Table 3. 



Extrapolation of Square End 


Stiffness 


end trays 
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TABLE 3. TRAY PANEL EFFECTIVE ORTHOTROPIC PROPERTIES 


Type 

t 

t 

t 

d 

d 

d 


u 

Side Trays 

X 

cm 

c^a 

cS' 

X 

cm 

cl 

c^ 

yx 

xy 

Debris 

.0493 

.0112 

.0211 

9.096 

9.627 

4.049 

-.116 

-.512 

Shallow 

.0500 

.0109 

.0142 

9.174 

9.716 

4.016 

-.118 

-. 538 

Nominal 

.0950 

.0282 

.0254 

9.136 

9.474 

-1.199 

-. 105 

-.351 

Deep 

.100 

.0323 

.0188 

9.098 

9. 558 

-25.077 

-.111 

-.344 

End Trays 
Rectangular 

Shallow 

. 104 

. 121 

.0254 

7.836 

7.859 

1.803 

-. 0742 

-.0641 

Nominal 

.113 

.129 

.0198 

7.922 

7.963 

.305 

-. 0461 

-.0404 

Deep 

.121 

. 136 

.0145 

7.996 

8.057 

-2.316 

-.0217 

-.0193 

Square 

Shallow 

.113 

. 113 

.0254 

7.849 

7.849 

1.803 

-. 0692 

-.0692 

Nominal 

.121 

. 121 

.0198 

7.943 

7.943 

.305 

-. 0433 

-.0433 

Deep 

.129 

. 129 

.0145 

8.026 

8.026 

-2.316 

-.0205 

-.0205 
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STATIC TEST - ANALYSIS EVALUATION OF ATTACHMENTS 


Static test data was used to evaluate the SPAR model longeron-to- 
intercostal joint condition and the diagonal end-fitting flexibility. The two 
types of static tests that were performed are shown in Figure 24. The SPAR 
static deflections are shown as a percent of the test deflections for both the 
bending and torsion tests. The model and tests were performed on the LDEF 
configurations without trays. The ball and pin joints refers to the longeron-to- 
intercostal joints previously discussed and shown in Figures 5, 12, and 13. From 
this information it was decided to perform first the dynamic analysis without 
trays, using the model that included the diagonal end flexibility. The pin joint 
model was considered to be the intent of the longeron-to-intercostal joint de- 
sign. However, it was observed that there was some out of plane motion 
because of the tolerance in the clevis joints that led to the consideration of 
modeling using ball joints. The true joint condition is someplace between 
these two. For the dynamic analysis without trays, it was decided to evaluate 
both the ball joint and pin joint conditions for comparison because the static 
test deflection correlation with analysis was not conclusive. 
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LDEF Configuration Without Trays 



Torsion 




Without End 

With End 

Without End 


Flexibility 

Flexibility 

Flexibility 


(%) 

(%) 

(%) 

Ball 

Joints 

77.7 

99.7 

94.4 

♦ 

Pin 

Joints 

71. 1 

89.6 

Percent 

75.6 

SPAR deflection 

= X 100 

test deflection 

* 

Longeron-to-Intercostal Joints 

** 


End Flexibility Refers to the DiagonaHo-Center Bhd Flat Plate Fittings 


With End** 
Flexibility 
(%) 

110.6 

84.8 • 


CJi 

05 


Figure 24. Analytical and Measured Static Deflections 



MASS PROPERTIES 


As previously discussed in the section on the SPAR LDEF model, the 
model was divided into nine substructure groups as shown in Table 1. For each 
SPAR group in Table 1, the consistent mass matrix is generated by SPAR and 
the weight of each group is printed out. Table 4 shows the weight of eight of the 
nine groups representing the weight of the SPAR model without trays before the 
diagonal end flexibility beam elements were added. All of the groups had E21 
beam elements and the lECM shelf also had E43 panel elements. The total 
SPAR weight at the group level was compared with the drawing weight and the 
lumped masses in SPAR were adjusted at the group level for better agreement 
with the drawing weight. This procedure results in an accurate overall weight 
distribution control. The total SPAR weight without trays is 3648 kg (8041 pounds). 
The keel pin and the damperMdghteare not included in this model or the test con- 
figuration. This model weight was 0. 2% lighter than the weight statement. The 
final lumped masses and grid point locations for the final configuration with trays 
are tabulated in Table 5. Most of these weights are the end fitting weights and 
miscellaneous hardware. Table 6 shows SPAR weight distribution for the final 
configuration with trays. The total SPAR tray weights, 5218. 1 kg (11504 pounds), 
exactly equals the test configuration and includes the weight of the trays, the base 
plates and the dummy weights that were bolted to the base plate. These total 
weights for each tray are tabulated in Figure 21. The total with tray SPAR model 
weight was 9057 kg (19967 pounds). 

The total weight for each tray was distributed uniformly over the 

surface of the effective panel area. The radial location of the tray and 

base plate center of gravity is based on SPAR weight outputs of the static 

tray models. The effective panel area of bays A, B, E and F was 
2 2 

1.4921m (2312.7 in ). The two center bays C and D had an effective panel 
2 2 

area of 1. 6197 m (2510.6 in ). The mass moments of inertia of each tray 
due to the radial offset of each tray center of gravity from the plane of the 
four comer longeron grid points was determined and one fourth was applied 
to each of the four comer grid points as shown in Appendix C. 
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TABLE 4 . SPAR GROUPING FOR WEIGHT BALANCE, 
CONFIGURATION WITHOUT TRAYS 


SPAR 

Group 

Title 

SPAR 

Lumped 

Total 

Drawing 


Number 


Weight 

Mass 

SPAR 

Weight 




(kg) 

(kg) 

(kg) 

(kg) 

E21 

1 

Longeron 

910.3 

189.1 

1099.4 

1082.6 

E21 

2 

Intercostals 

355.9 

0.0 

355.9 

366.7 

E21 

3 

Center Bhd. 

639.5 

278.7 

918.2 

945.6 

E21 

4 

Forward Bhd. 

188.8 

64.0 

252.8 

257.0 

E21 

5 

Aft Bhd. 

188.8 

64.0 

252.8 

257.0 

E21 

6 

Diagonals 

279.0 

154.7 

433.7 

433.8 

E21 

7 

Cross Beam 

94.2 

75.3 

169.5 

169.4 

E21 

8 

lECM Shelf 

16.6 

4.5 

26.3 

26.2 

E43 

1 

lECM Shelf 

5.2 





Sub Total 


2678.3 

830.3 

3508.6 

3508. 3 

Misc. 

Hardware (on Bhds. ) 


— 

55.2 

55.2 

55.2 

2 Trunnion Instl. (on Center Bhd. ) 


56.5 

56.5 

57.3 

Spindle Pin Instl. (Fwd. Bhd. ) 


19.8 

19.8 

19.8 

80 Clips (Fwd. and Aft Bhd. ) 



7.4 

7.4 

7.4 




2678.3 

969.2 

3647. 5 

3648. 0 


cn 

00 



TABLE 5. LUMPED MASS LOCATIONS, 

CONFIGURATION WITH TRAYS 


Grid 

Point 

Mass 

Grid 

Point 

Mass 


(kg) 


(kg) 

1 

26.08 

79 

30.37 

5 

26.08 

80 

30.16 

9 

32.99 

81 

22.01 

10 

15.32 

82 

34.61 

11 

.66 

83 

78.24 

12 

. 66 

84 

50.07 

13 

.66 

87 

50.07 

14 

.66 

88 

78.24 

15 

15.32 

99 

2.25 

16 

15. 32 

100 

2.25 

17 

. 66 

149 

15.32 

18 

.66 

150 

. 66 

19 

. 66 

151 

.66 

20 

.66 

152 

.66 

21 

15.32 

153 

.66 

22 

1.65 

154 

15.32 

23 

1.65 

155 

15.32 

24 

1.65 

156 

.66 

27 

1.65 

157 

.66 

28 

115.43 

158 

. 66 

29 

1.65 

159 

.66 

32 

1.65 

160 

15.32 

33 

1.65 

161 

1.65 

34 

1.65 

162 

1.65 

71 

34.61 

163 

1.65 

72 

22.01 

166 

1.65 

73 

30.16 

167 

1.65 

74 

30.37 

170 

1.65 

75 

22.01 

171 

1.65 

76 

34.61 

172 

1.65 

77 

34.61 

181 

133.24 

78 

22.01 





Total 

1102.00 




(2429. 5 lb) 



TABLE 6. SPAR WEIGHT, CONFIGURATION WITH TRAYS 


SPAR 

Lumped 

Total 

Weight 

Mass 

SPAR 

(kg) 

(kg) 

(kg) 


Longerons 

910.3 

189.1 • 

1099.4 

Intercostals 

355.9 

0.0 

355.9 

Center Bhd. 

639.5 

278.6 

918.1 

Fwd. Bhd. 

188.8 

64.0 

252.8 

Aft Bhd. 

188.8 

64.0 

252.8 

Diagonals 

315.8 

154.7 

470.5 

Cross Beam 

116.1 

75.3 

191.4 

lECM Shelf 

21.7 

4.5 

26.2 

Miscellaneous 


271.8 

271.8 

Total Without Trays 

2736.9 

1102.0 

2828.9 

Trays 

5218.1 

0.0 

5218. 1 

Total 

7955. 0 

1102.0 

9057. 0 
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COMPARISON OF FREQUENCIES AND MODE SHAPES 


Configuration Without Trays 

The comparison of some of the SPAR and test vibration modes is 
shown in Figure 25 for the LDEF without trays. The ball-joint and pin-joint 
analyses bracket the test data at the lower LDEF frequencies. In general, the 
longeron-to-intercostal ball joint boxmdary condition was closer to the test data 
so this boimdary condition was selected for the vibration analysis with trays. 
Analytical frequencies were higher than test for the second vertical and lateral 
bending modes. There were many intermediate modes not showry especially in 
the 16 to 30 Hertz range where the predominate response was diagonal motion 
often coupled with longeron lateral motion. 

The first SPAR lateral bending mode shape at 20.85 Hertz is compared 
with the test mode at 21.38 Hertz in Figure 26. The dashed lines represent 
the undeformed test shape. The straight solid longeron lines are those that 
had! no transducers in the test. By comparing the shape of individual longerons 
and diagonals, reasonable mode shape comparisons exist. The first SPAR 
vertical bending mode shape at 28.39 Hertz is compared with the mirror image 
of the test mode at 28. 52 Hertz in Figure 27. Again the individual longeron 
and diagonal mode shapes compare. Figure 28 shows the first SPAR torsion 
mode at 12. 94 Hertz compared with Figure 29 showing the test mode at 
14. 14 Hertz. 

Configuration With Trays 

The LDEF configuration with trays was evaluated using the longeron-to- 
intercostal ball joint SPAR model. Figure 30 shows the comparison of a few 
SPAR and test modes of LDEF with trays. The SPAR and test frequencies 
of the first three fundamental LDEF modes had very close correlation con- 
sidering the complexity of the variation in trays and the lor^eron-to-intercostal 
joint motions. 
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Figure 25. LDEF SPAR 



Cross Beam Rigid Body Yaw 
Cross Beam Longitudinal Bending 
First Torsion 
First Lateral Bending 

Predominantly Diagonal Motion 

First Vertical Bending 
Second Torsion 
Second Vertical Bending 
Second Lateral Bending 


and Test Vibration Modes Without Trays 









(Longeron to Intercostal Ball Joints) 


Frequency, Hz Mode Description 


SPAR 

Test 


5.72 

6.10 

Cross Beam Rigid Body Yaw 

9. 33 

10.10 

Cross Beam Longitudinal Bending 

19.06 

19.20 

1st Torsion 

22.75 

21.83 

1st Lateral Bending 

25.74 

23.90 

1st Vertical Bending 


Figure 30. SPAR and Test Modes With Trays 



REMARKS AND CONCLUSIONS 


The LDEF structure presented some unique features to challenge the 
finite element modeling techniques. 

1. The design of the longeron-to-intercostal clevis pin joint modified 
the LDEF bending stiffness and torsional stiffness. Static bending and tor- 
sional test data of the LDEF vehicle without trays installed was valuable in 
assessing this joint boundary condition. The ball-joint boundary condition was 
selected as more accurate than the pin-joint boundary condition based on this 
static test data and by comparing the vibration analysis of the two boundary con- 
ditions with vibration test data of the LDEF configuration without trays. 

2. The design of the diagonal end fitting that connected two diagonals to 
a common point on the center bulkhead was a tee shaped flat plate. This fitting 
appeared to be stiff enough so that it could be treated as a rigid link. However, 
there was a slight angle offset between the diagonal and the plane of the end fitting 
so that axial loads in the diagonal resulted in bending and twisting of this fitting. 
Again, the static bending and torsional test data was found to be valuable in deter- 
mining that the flexibility of this diagonal end fitting was required. 

3. The design of the tray-to-structure clamps shown in Figures 7 
and 8 allowed in plane relative sliding motion between the tray flange and 
the mounting structure. At the four side mid point clamps this relative 
sliding motion was approximated in the tray effective panel stiffiaess matrix 
by allowing degree of freedom 6 to be free when the tray finite element model 
was subjected to static shear and twisting deformations. This boundary change 
reduced the shear stiffness to 54 percent and twisting stiffiiess to 58 percent 

of the shear and twisting stiffnesses for the original boundary condition. 

4. Using the SPAR LDEF model and incorporating the modeling 
features described above resulted in good agreement between analytical and 
measured data. This was true for static deflections and vibrations modes. 
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APPENDIX A 


PANEL AND BEAM PROPERTIES 



TRAY PANEL STIFFNESS 

Effect of Boundary Conditions on Tray Panel Stiffness 

After examing the clamping conditions of the tray lip to the mounting 
structure it was decided to allow some clamp sliding due to shear loading. 

As shown in Figure 8, the clamp friction is the only restraint for relative 
in-plane motion, since there are no through bolts. The 4 clamps at the 
center of each side of the tray were allowed to slide parallel to the edge of 
the tray to relieve some of the shear loads being transferred at these clamp 
locations. No sliding was allowed at the 4 comer clamp locations. By 
freeing the degree of freedom 6 boundary condition at the ends of four rigid 
links as shown below in the right hand column, it was determined that the 
total shear transfer per side was reduced to approximately one half when 
the unity shear strain loading conditions were applied. The Debris tray 
stiffness matrix was derived with and without the relaxed boundary condition 
being used for the unit shear and unit twist cases. The numerical comparison 
is shown in Table 7. The shear and twisting stiffness terms were approxi- 
mately one half as stiff with the relaxed boundary condition. This relaxed 
boundary condition was used for the derivation of all of the effective panel 
stiffness terms with the values tabulated in Table 8. 

Original Boundary Relaxed Boundary 

(No Relative Tray Lip Motion) (With Relative Tray Lip Motion) 

(See Table 2) 

Grid 0 

Point No. 

1 0 

2 0 

3 0 

4 0 

5 0 

6 0 

7 0 

8 0 

SPAR BEAM SECTION PROPERTIES 

The SPAR beam configurations and section properties are tabulated in 
Table 9. The locations of these beam elements in the SPAR model are indicated 
in Figures 11, 14 and 19. 


0 

Free 

0 

0 

Free 

0 

Free 

Free 
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TABLE 7. FLANGE ROTATION EFFECT ON 

DEBRIS PANEL STIFFNESS MATRIX 

Relative Flange Motion Allowed ( 6^ = Free) 



Relative Flange Motion Not Allowed (0 = 0 ) 

z 



where the xinits are 



The only stiffness terms that change are the C.. , and 

oo ob bo bb 

terms. By allowing the flange motion the shear stiffness, , reduced from 
.1034 to .05594 or only 541 percent as stiff. The twisting stiffness, Cgg , re- 
duced to 58.1 percent and the coupling term, C^g , reduced to 55.8 percent as 
stiff. 
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TABLE 8. TRAY PANEL STIFFNESS TERMS 


(Side Flanges Free to Rotate) 


Cll 

Ci2 

C22 

C33 

Cl4 

Cl5 

C24 

MN/cm 

(MN/ cm) 

(MN/cm) 

(MN/cm) 

(MN) 

(MN) 

(MN) 


Side Trays 








Debris 

.3604 

-. 04186 

.08170 

.05594 

3.278 

-.3314 

-.2072 

Shallow 

.3690 

-.04361 

.08105 

. 03748 

3.386 

-.3616 

-.2943 

Nominal 

.6805 

-.07112 

. 20245 

.06772 

6.217 

-. 5654 

-. 3527 

Deep 

.7192 

-.07972 

.2319 

. 04963 

6.544 

-. 5712 

-. 3671 

End Rectangular 








Shallow 

.7234 

-.05368 

.8375 

. 06713 

5.669 

.03661 

. 00934 

Nominal 

.7788 

-.03590 

.8877 

.05277 

6.169 

.08776 

. 06107 

Deep 

.8341 

-.01813 

.9380 

.03841 

6.669 

.1389 

.1128 

End Comer 








Shallow 

.7805 

-. 05401 

.7805 

.06713 

4.448 

-.4261 

-.4261 

Nominal 

.8333 

-.03604 

.8333 

. 05277 

6.618 

-. 2862 

-.2862 

Deep 

.8861 

-.01816 

.8861 

. 03841 

7.113 

-. 1458 

-. 1458 


C25 

C36 . 

C44 

C45 

C55 

C66 



(MN) 

(MN) 

(Mm-N) 

(Mm-N) 

(Mm-N) 

(Mm-N) 


Side Trays 








Debris 

.7864 

-.2264 

.2564 

-.01664 

.07187 

. 009610 


Shallow 

.7873 

-.1505 

.2678 

-.02001 

. 07249 

.006602 


Nominal 

1.918 

.08118 

.3868 

-.01580 

.1655 

.003100 


Deep 

2.216 

1.2446 

.3480 

. 00264 

.1866 

. 02307 


End Rectangular 








Shallow 

6.582 

-. 1210 

.4365 

. 01427 

.5481 

. 004519 


Nominal 

7.069 

-.01606 

.4802 

. 01271 

.5922 

.009850 


Deep 

7. 556 

.08892 

.5232 

.01115 

.6363 

. 02106 


End Comer 








Shallow 

6.126 

-.1210 

.4923 

. 01427 

.4923 

. 004519 


Nominal 

6.618 

-.01606 

.5362 

. 01271 

.5362 

.009850 


Deep 

7.113 

. 08892 

.5797 

.01115 

.5797 

. 02106 



TABLE 9. SPAR SECTION PROPERTIES 


Beam Elements 

SPAR 

Beam 


Section 

Property 

SPAR 

Designation 

B1 

T1 

B2 

T2 

B3 

T3 

Area 


^2, 

J 

BAl 

WFL 

(cm) 

11.43 

(cm) 

.9652 

(cm) 

15.8496 

(cm) 

.5588 

(cm) 

(cm) 

(cm2) 

30.92 

(cm4) 

240.4 

(cm^) 
1746. 7 

(cm^) 

7.8 

BA 2 

TEE 

11.43 

.9652 

7.6708 

.5588 

- 

- 

15.32 

120.2 

79.4 

3.8 

BAS 

WFL 

17. 78 

1.4732 

83.1514 

.9144 

- 

- 

77. 56 

1381.8 

12619.7 

44.9 

BA4 

WFL 

11.43 

.9652 

22.1996 

.5588 

- 

- 

34. 47 

240.5 

3471. 2 

8.1 

BAS 

GIVN 

35. 56 

1.270 

25.40 

.635 

8.89 

1.016 

129. 03 

12444. 1 

17975.4 

8495. 3 

BA6 

WFL 

10.16 

.762 

13. 716 

.4572 

- 

- 

21.75 

133.3 

910.5 

3.4 

BAT 

GIVN 

26.518 

.3556 

15.24 

.4572 

3.7338 

.762 

42.78 

3788.6 

1940. 5 

2080.6 

BAS 

TUBE 

6.985 

7.620 

- 

- 

- 


29.14 

778. 3 

778.3 

1556.6 

BA9 

TUBE 

3.81 

5.08 

- 

- 

- 

- 

35. 47 

357.6 

357.6 

715.1 

BAIO 

GIVN 

38.10 

1. 5748 

12.70 

1.27 

22.86 

2.54 

116. 52 

10570.2 

877.8 

1216.6 

BAll 

GIVN 

5.08 

3.302 

8.636 

.762 

2.54 

.508 

32.52 

119.3 

426.5 

63.7 

BA12 

CRN 

4.445 

.6604 

6.2992 

.4318 

- 

- 

8.59 

80.3 

17.2 

1.0 

BA13 

ANG 

2.235 

.1600 

2.380 

.1600 

- 

- 

0. 74 

0.7 

0.2 

0.01 

BA14 

TUBE 

2.416 

2.54 

- 

- 

- 

- 

1.94 

6.0 

6.0 

11.9 

BA15 

TUBE 

.254 

6.9088 

- 

- 

- 

- 

149. 75 

1789.4 

1789.4 

3578.7 

BA16 

GIVN 

( Same as BA2 for Ball Joints) 

- 

- 

15.32 

0 

0 

0 

BA 18 

GIVN 

( Same as BA2 for Pin Joints) 

- 

- 

15. 32 

120.2 

0 

3.8 

BA19 

GIVN 

18.288 

2.54 

- 


- 

- 

46.45 

1294.6 

25.0 

99.9 

BA20 

GIVN 

17.78 

1.905 

- 

- 

- 

- 

33.87 

892.3 

10.2 

41.0 


M 


1 .) I 


T 


-3 
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JL 
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X 

Nl 

CQ 
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The bending stiffness for a single panel is: 

3 


D = 


Eh 
12 (1 - 


^ t = . 3556 cm (. 14 in) 

d = 15.24 cm (6. in) 

± L. 

t = . 3556 cm (. 14 in) 


T 


This bending stiffness was derived in Reference 2 from the integration 


h/2 

M = ^ z d z 
-h/2 

For a two plate configuration the integration is 

-d/2 d/2 + 1 

M = \ cyzdz+ \ azdz 

J X J X 

-d/2 -t d/2 

resulting in the bending stiffness 
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APPENDIX B 


MESH SIZE CRITERION 



MESH SIZE CRITERION 


Maximum finite element grid intervals or mesh sizes were determined 
by use of a criterion stating that a beam or plate element having simply sup- 
ported boundary conditions should have a natural frequency greater than 50 
Hertz. If the natural frequency was less than 50 Hertz, intermediate grid 
points were added until the criterion was met. The value of 50 Hertz was 
selected because the region of LDEF frequencies of interest were those up 
to 50 Hertz . As shown in Figure 31, the only place that intermediate grid 
points were required was on the diagonals and lateral members of the cross 
beam. 
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Required: Frequency > 50 Hz 


Boundary Condition: Simply Supported 

B eam Elements 

Max. Grid ^ 

Interval 

S 

Failed Criteria 
( < 50 Hz ) 

(Intermediate grid points added) 

Longerons Diagonals 

Intercostals Cross beam - lateral elements 

Bulkheads 
lECM Shelf 

Cross beam - vertical center 
element 

Plate Elements 


r 

tSr 

Met Criteria 
( > 50 Hz) 


ss 


ss 

Max. Mesh 
Size 

ss 


ss 


Met Criteria 
( > 50 Hz) 

All plate elements met the criteria 


Figure 31. Mesh Size Criterion 
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APPENDIX C 


TRAY MASS MOMENTS OF INERTIA 



TRAY MASS MOMENTS OF INERTIA 

The mass moment of inertia of a given tray type (I^) was obtained from 
the SPAR tray substnicture group weights and the radial offset distance (z ) 

L 

to the tray mass center from the plane of the comer longeron grid points. These 

values are tabulated in Table 10. 

E wz 


_2 

I^ = z S w 

The average ballast height for each tray and the total ballast weight (w^) for each 
tray was used to determine the ballast offset distance (Zj^) and the ballast mass 
moment of inertia (Ij^). These values are tabulated in Table 11. 

2 

’b = S. *b 


The total tray radial offset mass moment of inertia for each tray is 


I 



One fourth of this value was then applied at each of the four comer grid points 
and rotated from the tray coordinate system to the LDEF cylindrical coordinate 
system. 

r z 



X, y, z Tray coordinates 
r, 0, z’ LDEF coordinates 


I 


0 


Mass Moments of Inertia in LDEF Coordinates 
2 2 

= I cos (p + l sin <0-2 1 sincocoscp 
y z ^ yz ^ ^ 
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I 

r 


2 2 

I sin 0 + 1 cos 0 + 21 sin 0 cos 0 
y ^ z ^ yz ^ ^ 

'er 

= 

I cos 20 + i (I -I) sin 20 
yz ^ 2 'y z' ^ 

For I 

y 

= 

1/4 ; I = I = 0. 0 
z yz 

'e 

= 

(l/4 ) cos^ 0 

'r 


( 1 / 4) sin^ 0 

^xe 


(l/ 8 ) sin ^0 

1 1 

Z 

= 

1/4 

ip 


± 15° 

2 

cos (p 


.93301 

. 2 

sin (p 

= 

.06699 


TABLE 10 » SUMMARY OF SIDE TRAY OFFSETS 
AND MASS MOMENTS OF INERTIA 

h 



Tray and 

Tray 


Tray 

Base Plate 

Mass 


Type 

Mass 

Offset 

It 


kg (lb) 

cm (in) 

ni^kg (in-lb- s 

Debris 

24. 05 (53. 023) 

2.78 (1.096) 

.01863 (.1648) 

Shallow 

26.09 (57.519) 

2.35 (.925) 

. 01443 (. 1276) 

Nominal 

29. 00 (63. 937) 

-3.53 (-1.390) 

.03614 (.3196) 

Deep 

32.24 (71.083) 

-9.01 (-3.547) 

.2619 (2.316) 

Deep* 

48.14 (106.137) 

-40.4 (-15.9) 

7.847 (69.40) 


Special deep tray at location D-12 



TABLE 11. SIDE TRAY BALLAST MASS MOMENTS OF INERTIA 








Zb 


Tray 

Ballast 

Ballast 


Tray 

Ballast 

Ballast 


Location 

Weight 

Offset 


Location 

Weight 

Offset 



(kg) 

(cm) 

(m^-kg) 


(kg) 

(cm) 

(m^-kg) 

A-1 

0 

- 

0 

D-1 

73.94 

4.902 

.1777 

A -2 

70.31 

-3. 037 

.06486 

D-2 

0 

- 

0 

A-3 

51.71 

4.318 

. 09643 

D-3 

0 

- 

0 

A -4 

68.04 

-3. 037 

.06278 

D-4 

73.94 

4. 902 

.1777 

A-5 

0 

- 

0 

D-5 

0 

- 

0 

A-6 

0 

- 

0 

D-6 

70.31 

-3. 037 

. 06486 

A-7 

16.33 

2.945 

. 01416 

D-7 

0 

- 

0 

A-8 

0 

- 

0 

D— 8 

74.84 

4.903 

.1799 

A-9 

70.31 

-3. 037 

.06486 

D-9 

0 

- 

0 

A-10 

70.31 

-3. 037 

. 06486 

D-10 

27.22 

3.326 

. 03012 

A-11 

0 

- 

0 

D-11 

73.94 

4. 902 

.1777 

A-12 

38.56 

3.835 

.05672 

D-12 

131.09 

-49. 530 

32. 159 

B-1 

0 

- 

0 

E-1 

16.33 

2.945 

. 01416 

B-2 

68.04 

-10.659 

.7730 

E-2 

68.04 

-10.659 

.7730 

B-3 

70.31 

-3. 037 

.06486 

E-3 

65.32 

4.901 

.1569 

B-4 

8.16 

3.254 

.00864 

E-4 

70.31 

-3. 037 

. 06486 

B-5 

68.04 

-10.659 

.7730 

E-5 

68.04 

-10.659 

.7730 

B-6 

13.15 

3.125 

. 01284 

E-6 

28.12 

3.327 

. 03112 

B-7 

68.04 

-10.659 

.7730 

E-7 

68.04 

-10.659 

.7730 

B-8 

8.16 

3.254 

.00864 

E-8 

70.31 

-3. 037 

. 06486 

B-9 

70.31 

4. 903 

.1690 

E-9 

70.31 

-3. 037 

. 06486 

B-10 

31.75 

-12. 182 

.4712 

E-10 

68.04 

-10.659 

.7730 

B-11 

0 

- 

0 

E-11 

0 

- 

0 

B-12 

8.16 

3.254 

.00864 

E-12 

10.43 

-12.616 

.1660 

C-1 

8.16 

3.254 

.00864 

F-1 

16.33 

2.945 

. 01416 

C-2 

0 

- 

0 

F-2 

11.79 

2.996 

. 01058 

C-3 

26.76 

3.250 

.02827 

F-3 

37.19 

-4. 181 

. 06500 

C-4 

27.22 

3.326 

. 03012 

F-4 

0 

- 

0 

C-5 

42.18 

3.887 

.06373 

F-5 

11.79 

2.996 

. 01058 

C-6 

70.31 

-3.037 

. 06486 

F-6 

10.43 

-4. 946 

. 02551 

C-7 

0 

- 

0 

F-7 

0 

- 

0 

C-8 

73. 94 

4. 902 

.1777 

F-8 

44. 00 

3.886 

. 06644 

C-9 

21.77 

3.125 

. 02126 

F-9 

37. 19 

-4. 181 

. 06500 

C-10 

8.16 ■ 

3.254 

.00864 

F-10 

0 

1 

0 

C-11 

73.94 

4. 902 

.1777 

F-11 

0 

- 

0 

C-12 

83.57 

3.632 

. 04428 

F-12 

44.91 

.-4.- 054 

.07382 


APPENDIX D 


SPAR EIGENVALUE EXTRACTION 


EXPERIENCE 



SPAR EIGENVALUE EXTRACTION EXPERIENCE 


Vibration inodes are computed in SPAR by means of a sparse matrix 
eigensolver as described in reference 1. The procedure implements an 
iterative process consisting of a Stodola (matrix iteration) procedure followed 
by a Rayleigh- Ritz procedure, followed by a second Stodola procedure, etc. , 
resulting in successively refined approximations of a reduced set of m 
eigenvalues. Since K need not be positive definite, this method can be used 
to solve "shifted" vibration eigenproblems. This "shifted" procedure is 
required for a free-free system to make the stiffness matrix non-singular. 

The procedure used in SPAR is to alter the original linear vibration problem 

r Mx - Kx = 0 (1) 

by adding and subtracting cMx from the left side, 

(r - c) Mx - (K - cM) x = 0 (2) 

Equation 2 is of the same form as equation 1, except that (r - c) has re- 
placed r, and (K - cM) has replaced K. The value c is the shift value. 

SPAR inverts the shifted stiffness matrix (K - cM) and computes the 
eigenvalues in the neighborhood of c. When this inversion is being made, the 
number of negative terms encountered in factoring K - cM is equal to the 
number of roots (eigenvalues) below the shift value (c) in equation (2). 

Using this SPAR procedure for the LDEF model without trays 
required eleven computer runs to obtain all of the flexible modes up to 
50 Hertz as illustrated in Figure 32. The inputs shown for each computer 
run are the upper and lower frequency bounds (VI and V2) and the shift 
location indicated by the triangle. The outputs shown are the number of 
converged eigenvalues and the frequency range of these converged eigen- 
values. The first run is shown on the bottom right of the figure with a shift 
at 50 Hertz and the last run is shown at the top right. From the first run. 
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which was used to determine how many eigenvalues existed below 50 Hertz, 

93 flexible modes were indicated, and 15 converged modes were found. The 
next four runs directly above and to the left of the first run in Figure 32 were 
run simultaneously. From these runs, thirty-four additional converged modes 
were identified, leaving forty-four that were found in the remaining runs. 

At the lowest frequency range, only two modes converged. This 
is partially due to the fact that there were only two modes in the immediate 
vicinity of the shift frequency while the other runs had many closely spaced 
modes near the shift location. However, there were 7 modes between the 
frequency bounds of VI and V2 in this low frequency run so that more 
converged modes were ejqjected. More problems were encoimtered in 
obtaining converged modes at the lower frequencies when the lower boimd 
VI was near the rigid body (zero frequency) modes. 

In all of the runs, the value of INIT (number of initially chosen vectors) 
was set equal to 15. The value of NDYN (number of iterations) was increased 
from the default value of 8 to values of 13 to 16. With NDYN equal to the 
default value, the niunber of converged modes per run decreased considerably 
in many runs, and this procedure was not considered cost effective for the 
LDEF model with 1416 degrees of freedom. By examining the runs in 
Figure 32, the number of converged modes after 8 iterations (NDYN = 8) was 
tabulated and compared with the converged modes from the computer runs. 

The results are shown as follows; 


Shift 

Frequency 

(Hz) 

NDYN 

in 

Run 

No. of 
Converged 
Modes from Run 

No. of Converged 
Modes if 

NDYN = 8 (default) 

7.96 

16 

2 

2 

13.96 

16 

12 

12 

18.15 

15 

6 

6 

20.00 

13 

12 

9 

22.51 

15 

12 

6 

29.00 

15 

10 

4 

35. 59 

16 

13 

13 

38.98 

16 

10 

5 

41.20 

15 

4 

4 

43.59 

16 

11 

6 

50.00 

16 

15 

9 
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The increased number of converged modes obtained with an increased 
number of iterations may be due to the large number of degrees of freedom 
and the numerous modes with nearly identical frequencies. 
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Legend 
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Figure 32. SPAR Eigensolutions for Free- Free Modes of LDEF 
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